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We consider tensor-multiscalar representations for several types of modified gravity ac- 
tions. The first example is the theory with the action representing an arbitrary smooth 
function of the scalar curvature R and OR, the integrand of the Gauss-Bonnet term and 
the square of the Weyl tensor. We present a simple procedure leading to an equivalent the- 
ory of a space-time metric and four auxiliary scalars and specially discuss calibration of a 
cosmological constant and the condition of the existence of dS'-like solutions in the case of 
empty universe. The condition for obtaining a smaller number of independent scalar fields 
is derived. The second example is the Eddington-like gravity action. In this case we show, 
in particular, the equivalence of the theory to GR with the cosmological constant term, with 
or without use of the first-order formalism, and also discuss some possible generalizations. 
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I. INTRODUCTION 



Recently, there was a considerable interest in the f{R) theories described by the action 

= j d^x^gf{R), (1) 
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where f{R) is some difFerentiable function, see [1] for recent reviews. It is well-known that, under 
the condition f"{R) ^ 0, the theory is dynamically equivalent to scalar-tensor theory of gravity 
with the potential depending on the form of the function f{R). Our purpose is to discuss this 
equivalence in a slightly different framework. As an application of our method we will be able to 
generalize the equivalence theorem to the more general case when the action depends on a function 
of many variables, /(X*), with being, e.g., R, OR, Gauss-Bonnet integrand E and/or other 
quantities. Let us note that quantum corrections to GR (coming from the semiclassical approach 
to quantum gravity or from the string theory) can be modelled by such a function to some extent. 
Another advantage of the method which we present here is that it can be used, also, for other 
theories, e.g. for Eddington-like models. Last, but not least, our method is a bit more explicit and 
simple than the previously known ones, e.g., [2-8]). All considerations will be presented for the 
D = A case in order to make them more explicit, but they can be more or less straightforwardly 
generalized for an arbitrary D ^ 2 case. 

In the present paper we will mainly follow Lagrangian approach, but it is worthwhile to mention 
that there are also some recent papers treating f{R) [9] and even more general f{R^vap) [10] 
theories within the canonical formalism. Earlier, the discussion of stability issues for a while class 
of higher derivative theories has been presented in [11] at both classical level and also for the 
Wheeler-DeWitt equation. In general, the equivalence of two theories at the classical level does 
not imply their equivalence at the quantum level. The quantum treatment of higher derivative and 
scalar-tensor theories is not a trivial issue, hence here wc present a short discussion of it only. 

The paper is organized as follows. In Sect. 2 we consider the simplest case of the theory (1) and 
describe a simple way of mapping it into a metric-scalar (scalar-tensor) model. The content of this 
section is mainly not original, we just give a bit more simple form of the known transformations. 
One of relatively new aspects of our consideration is the procedure to fix the cosmological constant 
term in the metric-scalar representation of the theory. We also check this procedure by using the 
dS-like exponential solution. In Sect. 3 we address gravity theories of a rather generic form, in 
which the action contains an arbitrary function of various scalar, curvature-dependent invariants, 
such as the scalar curvature R, the Gauss-Bonnet term, the square of the Weyl tensor and others. 
We develop a systematic approach to map such theories into the metric-multiscalar models. One 
of the new elements of our consideration is that it includes the case when the initial theory has 
such Hamiltonian constraints that the number of independent auxiliary scalars is smaller than the 
number of initial curvature-dependent invariants. The condition for this to occur is obtained. Sect. 
4 is devoted to the formulation of general conditions for the existence of exponential solutions in 
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the theories investigated in Sect. 3. In Sect. 5 we consider the theory which is based on the string 
low-energy effective action of gravity, up to the third order in curvature invariants. In this case 
one cannot construct an equivalent scalar-tensor representation at the level of action or general 
equations of motion for the metric field. However, this problem can be perfectly addressed for 
the much more restricted case of a homogeneous and isotropic cosmological solution. In Sect. 6 
we apply our method of constructing equivalent theories to the wide class of the Eddington-like 
gravity theories. It is shown how to construct dual theories for such gravity theories, including 
both the second order and Palatini formalisms. In the latter case the auxiliary tensor field can 
be interpreted as a space-time metric. In Sect. 7 we discuss the equivalence between different 
representations of the theory (1) and its generalizations at the quantum level. Finally, in Sect. 8 
we draw our conclusions. 



II. EQUIVALENCE OF f{R) AND METRIC-SCALAR THEORY 

Let us start from a simple pedagogical example of the theory (1) and find its metric-scalar 
(scalar-tensor) dual. Consider the theory described by the action 

Si = J d^x^{i/jR-Vm. (2) 

The theory (2) describes a dynamical scalar ip, despite there is no kinetic term for i/j in the action. 
One can establish the relation between the theories (1) and (2). The equation of motion which 
follows from the variation of -0 in (2) has the form 

iJ = F'W = ^. (3) 

After solving (3) with respect to and substituting this solution ip = ip{R) back into (2), we obtain 
the action (1) with 

i;{R) ■ R - V{^iR)) = f{R) . (4) 

This means that the equivalence of the two actions is dynamical, i.e. it holds on extremal curves 
of the field tp. Later on we shall confirm the validity of this procedure through the equations of 
motion for both metric and i.e. in a way similar to the one of [5]. 

Our next step will be to find the relation between the functions V{'il>) and f{R)- Taking the 
derivative d/dR of Eq. (4), we arrive at the relation 

+ Ri^'iR) - V'iiP) ^'{R) = f\R) . (5) 
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In this formula we assume that = ip{R) and R = R{i^)- Using (3), the equation (5) immediately 
reduces to the very simple relation 

^ = f'{R), (6) 

indicating that the function R = V'{tp) is nothing else but the inverse to the function ip = f'{R)- 
Finally, we arrive at the following receipt for deriving the potential V{ip) for a given f{R)- 

a) Calculate = f'{R) and invert it, obtaining R = V'{iIj). Note that the possibility of such 
inversion requires f"{R) 7^ 0. 

b) Integrate over ip: 

v{i^) = no+ r R{i^)d^. (7) 
Jo 

One has to note that an arbitrary integration constant Qq in (7) exactly corresponds to the 
constant /o = f{R = 0) component of the integrand of Eq. (1), which is indeed lost when we 
take the derivative f'{R)- Furthermore, in order to fix the constant Qq, one can use the following 
simple consideration. By using (4) we arrive at 

V{i)) = Rip- f{R), where ^ = f{R). (8) 

Remember that when placed into the covariant action, can not be regarded as an irrelevant 
constant, because it is multiplied by the metric-dependent factor ^/^. As far as (7) should be 
equal to (8), one can then fix Qq. Later one we will additionally check the validity of this procedure 
for a cosmological dS'-like solution. 

The prescription given above enables one, in principle, to find the potential function V{iJj) for 
a given f{R). Let us check the results of this simple procedure at the level of the equations of 
motion. Taking variation of the equation (1) with respect to the metric, we obtain 

f'{R^,-^Rg^,) + ^g^,(Rf'-f)-V^V.f' + g^.nf = 0. (9) 

Performing the same operation for (2), we arrive at 

V'(i?;,.-^i25M-) = ~9^^uVii^) + - g^^n^ = 0. (10) 



One can verify that the equivalence between (9) and (10) holds if the relation 

f'R-f ^ YM. 

is satisfied. It is easy to check that the solution of this equation has the form (6). 



(11) 
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Consider some particular example for the procedure described above. The simplest case leading 
to the linear equations ijj = f'{R) is 

f{R) = n-K'^R + ^R\ (12) 

Using our previous results, one can easily arrive at 

M) = f'iR) = -K^ + aR =^ y'(V) = R = t±J^_ (13) 

a 

Integrating (13) we get 

Finally, in order to fix the integration constant, one has to put ^p = —k^ + aR back into (14) and 
compare it to (12). This procedure gives us 

n = no- — . (15) 

One can perform a simple verification of the described procedure for fixing J7. For this end, we 
will now derive the dS-\ike solution for both theories (1) and (2) in case of (12). The metric of our 
interest is 



a{t) = exp{cr(t)} and afterwards we will set a{t) = Hot. It is easy to obtain the equation for a for 
the theory (2), 

1 SS 

— ^ = -6 e-2'^(2i/' k + ip" + 2a" il) + 2il;'a' + 2a'^ip) - 4V{ip) = . (16) 



5a 

Here the prime stands for the derivative with respect to the conformal time, e.g., 

, da , s. da 
a' = — = ait) — , 
df] ^ ' dt' 

while the derivative with respect to the physical time t is denoted as a dot. It terms of the physical 
time and adopting A; = 0, we obtain the relation 

-12i/^V-9i?o^-3V' = 2V . (17) 

Using ip = —K? + aR and taking into account R = —12 Hq for the FRW metric, we get 

V = GH^K^ -6H^aR (18) 
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and finally 

n = QHIk^ . (19) 

On the other hand, starting from 

Sf = J (fx^f{R) = j d'^x^ n-K'^R+^R^ (20) 

we arrive at the equation 

^ ^ = _ 6Ac2e-2<^(2a'2 + 2a") + ISae'"''' {2a"" - 12a' V") = . (21) 

It is straightforward to check that the solution a = Hot = — ln(iJo|?7|) corresponds, again, to the 
relation (19). 

Finally, for the sake of completeness, let us address the possibility of using conformal trans- 
formation to deal with the metric-scalar theory. It is well known fact that the theory (2) can be 
easily mapped into another one with the standard form of the scalar kinetic term. We will give 
the corresponding treatment here just for completeness and refer the reader to the review [12] for 
further details and (numerous, indeed) references. 

Let us start from the conformal transformation 

9^,u 9^. = 5m- • e'^^^) (22) 
in the action (2). Simple calculation yields the following result 

3^9^.^", i^]= J d^xV^ e^- [R - 6{Vaf - GOa] - e"'' F(V')} , (23) 

where (Vcr)^ = g^^d^adya. Let us choose a such that il^e^'^ = —k^. Then the first term if) e^'^ R 
coincides with the Einstein-Hilbert term, also the third term —6^ e^^Oa becomes a total derivative 
which does not affect the equations of motion. In order to provide the standard form of the kinetic 
term, we take 

ip = 2V^Ka, then ^ = -/^^ exp { p-}. (24) 

The output looks like 

Smin[9,.u, = J d^xV^ {-^R+1 gf^-'d^ifd,^ - U{^)} , (25) 
where the two potentials are related as 



K 
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U{ip) = e^'^Vii;) = ^V{^). (26) 
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The formulas (22), (24) and (26) are nothing but the change of variables in the action (2) which 
transform it into the action (25). Therefore we do not need to check the equivalence between 
two actions by other means, e.g. examining the equations of motion. Finally, let us note that 
the multiscalar case can be, in principle, also treated by the conformal transformation, but this 
transformation is not so easy as in the one-scalar case [13]. 



III. MODIFIED GRAVITY THEORY OF A MORE GENERAL FORM 



We have presented a useful and simple prescription of mapping theories (1) into theories (2) 
at the classical level. This method can be generalized to the gravitational actions which are more 
general than (1). However, as we will see in what follows, in this case one needs more scalar fields. 
Some similar results has been recently obtained in [14] and [15], but the derivation there looks 
rather different. 

The method may be especially useful for working out the dS-type solutions and, therefore, 
is applicable for testing various models of modified gravity, including the ones corresponding to 
quantum corrections. Consider the following gravitational action 



Sgen 



j d^x^f{R,UR,C\E), (27) 



where E = E — "^U R^ E = R^^y^p ~ ^^ai3 + the Gauss-Bonnet topological term (Euler 

density) and = R^^^^ - 2Rl^ + (1/3) R^ is the square of the Weyl tensor. In view of the 
cosmological applications, it proves more useful to consider E rather than E. 
We start by introducing a generalization of the action (1), 



J (fx^/^J{Xi) where Xi = R, UR, E, (28) 

in case of the action (27), but the number of invariants can be easily extended. For this end we 
set i = 1,...,N. Consider the dual action 

S2 = J d^x^[Xi^p' -V{^p')\, (29) 

where repeated indices imply summation, as usual. Let us follow the same scheme which we applied 
in the previous section. The equations for the ■0* have the form 

dV 

We put them into (29), demanding equivalence to the action (28), 

S2 = J d^arV=?{^^|^-F(VO}= J d'^xV^fiXi). (31) 
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Assuming 

fiXi) = Xi^p' - V{^P') (32) 
and taking partial derivatives with respect to in (32), we arrive at 

where we used Eq. (30). The formulas (30) and (33) show that we always have 

;.,.f?|M = (34) 

dil)i dXk 

After all, the prescription for deriving V{il^^) is very similar to the one described in the previous 
section and looks as follows: 

a) calculate ip^ = ^q^J ; 

b) solve these equations and find Xk{'^) = 

c) integrate the last relations and find F(^fe) up to the additive constant; 

d) fix this constant by the requirement that the actions coincide in the corresponding limit 
(typically zero curvature). 

One can note that this procedure can be applied also to the non-Riemannian generalizations of 

GR, including the theory of gravity with torsion. 

Let us consider an example of how the equivalent metric-scalar theory can be achieved. We 
start from the action (27) with the function 

fiR,E) = F{R)-E, (35) 

where F{R) is an arbitrary function of scalar curvature, which will be fixed later on. The equivalent 
action is 

S2 = J d^x^ {^i? + xE- F(^, x)} . (36) 
Let us follow the prescription described above. The equations 

IL^X = F(R), g^i,^EF-(R) (37) 
can be solved with respect to the two scalar fields, 
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On the other hand, we have inverse functions 



Then 



^(V',X) = J 9ix)dip + gi{x) = 5i(x)+V'5(x) (40) 



and 



If we compare the two forms of the potential function (40) and (41), it becomes clear that g{x) 
satisfies the functional equation 

and, moreover, 

g2{i^)-Ci; = gi{x) = Ci = const. (43) 

It does not look possible to advance further, so let us take a more concrete form of F{R). 

Consider first a very simple case f{R, E) = — + RE. Making the same steps as in the general 
case, we obtain 

F{R) = R^^ = R^ g{x) = x, F'{R) = 1. 
Using (40) and (41), one obtain 

Vii', X) = 9i{x)+i^X = y dx + 92ii') = ipx + 92ii')- 

It is easy to see that in this case gi{x) = 32 (V') = C, so we get 

V{^,X) = ^X + C. (44) 

Finally, inserting the relations x = -R and i/j = E into (44), one can easily verify that C = CI. This 
result can be also checked by inspecting exponential solutions in the two cases. We avoid to bother 
the reader with the details of this verification, but just note that its output is positive. 
Consider a bit more complicated example when 

R \ 



f{R,E) = n-K'^R + P E In (1 + —^ (45) 
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where fi is some constant and Rq\s & reference value for the scalar curvature. The expression (45) 
can be seen as a part of the renormalization group corrected vacuum action, where the renormal- 
ization group parameter /x^ is associated to the scalar curvature (see, e.g., [16] for further details 
and references). 

Following the footsteps of the previous examples, we derive 

^ «.;^(e./^-l) (46) 

and 

'-^-'^-Wk-^' ^ (47, 

At the next stage we find 

^^g^'^^ = R ^ mx) = 5i(x) + i?o^e^/'3-i2o^. 

= E^ F(V', x) = 92W + Ro^ e^/^ + 1^ Ro e^'^ . 
Using these two expressions it is easy to figure out that 

gi{x) = K^Roe^/^ + C and 52(V') = -i^oV' + C- 
Finally, we arrive at the potential 

F(V', x) = Roi^ (e^'^ - 1) + Ro e^/^ -VL- k'^Rq , (48) 

where the constant C = — O — k^Rq has been fixed following the same method which we used in 
the previous cases. 

One can consider more complicated expression for the covariant Lagrangian, 

f{R,E,C^) = n-K^R + /3Eln (^1 + + /5i In (l + • (49) 

At this point one can make an important observation. Prom the first sight, the equivalent La- 
grangian for this case should have three auxiliary fields, because there are three structures R, E 
and C^. At the same time, the problem of reducing the theory (49) is essentially equivalent to the 
one of the theory (45), with the /3 E traded by the combination /3 £J + /3i C^. Obviously, in this case 
we need only two auxiliary fields and not three of them. In other words, in this case the number 
of necessary auxiliary fields is smaller than the one which could be thought by just counting the 
number of the structures Xi in the starting action. This example shows that it would be interesting 
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to have a general criteria for establishing an exact number of necessary auxiliary fields for a given 
initial function f{Xi). 

The problem of our interest is very close to the one which is typical for the transition from 
Lagrange to Hamiltonian formalism in the theory with constraints [47] (see, e.g., well-known books 
[18, 19] for introduction purposes). Indeed, it is analogous to the passage from the Lagrangian 
description, with no explicit dependence on the coordinates, to the Hamiltonian one. In this case 
f{X^) and X^tpi — V{ijji) play the roles of the Lagrangian and the Hamiltonian respectively, where 
the "momenta" are defined by 

(50) 

Finally, the quantities of {X^,X'^,X^} play the roles of "velocities". For the specific case of eq. 
(49), these equations have the form 

^^ = -^+ R, + R ' (51) 
^-2 = /31n + (52) 



-^0 / 

= A In ( 1 + # V (53) 



^0 / 

from what we directly infer the presence of a constraint, which is a dependence relation between 
the V'i given by 

^ = M2 - = PiX -PX1 = 0. (54) 

Since the constraint ^(V'i) comes directly from the definition of the "momenta", it is classified 
as a primary constraint [18, 19]. We note that the equation ^(V'i) = defines a surface in the 
space {'01, V'2,'03}- It should be stressed that this constraint only represents a restriction in the 
"momenta" space, while in the "velocity" space {X^,X^,X^} the constraint does not lead to any 
restrictions, since if the Vi are written as functions of the X'^, the constraint ^(V'i) becomes the 
function 0(0j(X'^)), which is identically null. 

If this were a typical Hamiltonian problem, one would evaluate the evolution of the primary 
constraints (e.g., the single one in the example considered above) in search for further constraints. 
However our present problem is simpler, because it does not involve any dynamics. Hence only the 
constraints with no relation to dynamics are relevant here, which are the primary ones. 

Before proceeding towards the determination of the potential V, we remark here on the relation 
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between the Hessian matrix 



dX^dXi 

and the presence of constraints. Firstly, if the Hessian is non-singular, the inverse function theorem 
guarantees that (at least locally) one can use the definition (50) to express the as a function 
of the and thus no constraint is expected. However the Hessian can turn out to be degenerate. 
For instance, in the case of (49), this 3x3 matrix is a singular matrix of rank two (for R ^ —Ro), 
namely 



^-E/3 + C^/3i (5 

/3 

V /3i oy 



(55) 



Therefore, it has a single independent zero-mode (i.e., an eigenvector whose corresponding eigen- 
value is zero). Indeed, 

iy=(Q f}^ -fi) (56) 

can be promptly identified as the single linearly independent zero-mode of the Hessian (55). Let 
us note that we choose to work with the zero-modes that multiply the Hessian matrix by the left. 

To conclude this introductory part, we note that each independent zero-mode generates an 
independent constraint. In particular, by multiplying the zero-mode u on both sides of the definition 
(50), one finds the same constraint (54). Afterwards we will show that each independent constraint 
leads to an independent zero-mode of the Hessian matrix. We note that this simple relation 
between zero-modes and constraints does not have a counterpart in general Hamiltonian problems 
with constraints, in particular since the corresponding zero-modes may depend on "coordinates" 
there (in the present problem, we are considering the analogous Hamiltonian problem in which the 
Lagrangian only depend on the "velocities" X*). 

In the presence of constraints, one cannot use the first relation of (34) to find V, since this 
relation is not valid in the presence of constraints. Namely, consider the variation of V in the 
constraint surface (i.e., in the surface ^ = 0), 

5V = Sii^iX' -f) = X'di^i + (^Vi - ^) SX' = X'Si/ji . (57) 

In the last step above, we used the definition of the momenta (50). The previous equation shows 
that V can be written as a function of tp alone, even if constraints are present. Thus, using the 
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last equality, 



(^-X-|**=0. (58) 



Now, using the Theorem 1.2 of [19], we find the extension of the first relation of (34) to the 
constrained case. 

In this formula are Lagrange multipliers and c^^, with m = 1,2, ...,M, are all the constraints 
of the model under consideration. In the particular case of (49) we have a single constraint, i.e. 
M = 1. The introduction of these multipliers is necessary since the relation between the and 
the has to be extended in order to become invertible. 

Now we are in a position to discuss the method of constructing potential V{tpi) in the case of a 
theory with constraints. The integration method which was employed previously can be extended 
to this case. For instance, in the example of f{R, E, C^) theory (49) one can solve the definition 
of V'i and arrive at 

R = Ro ie^'^ - l] , (60) 



E = ^e^/0^^ + ^^)-^c\ (61) 
Hence, from (59) and (54) we find 

Ro{e^^^-l)^'^, (62) 
^e^/^{^ + n^)-^C^ = '^ + XP^, (63) 
C' = ^-XP (64) 

The first equation can be straightforwardly integrated, but the second cannot, since we do not 
know how to express both and A as functions of the tpi. Nevertheless, the Lagrange multiplier 
is still free, thus we can set it in such a way that the term disappear, namely A = —C'^/p. 
Consequently, the third equation (64) becomes 

After that the expression for V can be integrated, and V{ijji) can be found using procedures similar 
to that used in the unconstrained case. One should note that the number of auxiliary fields in 
this procedure is smaller than the number N of the Xi structures in the initial f{Xi) theory. For 



14 



example, in the (49) case, albeit we started from the f{Xi) which depends on three independent 
quantities, the corresponding potential V only depends on the two independent scalar fields. 

The form of the potential V can have some impact on the physical consequences of a given 
theory f{Xi). Hence, before concluding this example, we present y in a more general form. Let 

us start, as usual, from the simple example. In Eq. (60) we have made a choice of selecting the 
particular expression of i? as a function of x and xi- At the same time one can express R in a 
more general form, 

R = [e^//3 _ 1 + ^(eXiM _ 1)] , (65) 

where ^ is an arbitrary real number different from —1. The choice in (60) corresponds to ^ = 0. 
Prom the more general version (65) and Eq. (51) one can easily obtain 

PE + = i?o^^^ {e^/^ + ee^^/^^) . (66) 

Now, instead of attempting a direct integration of V (as in the non-constrained case) , we express 
it on the constraint surface ^ = 0, as 



(ex//3_i) + ^(exi//3i_i) 



^ - k'Ro- '-j^^ + Ex + C'xi 



(67) 



= + ^')- - ^ (68) 

= Ro{^ + Ac^) (e^/f^ (69) 

As one should expect, on the constraint surface V is independent on the value of ^. Nevertheless, 
a choice of ^ = 1 or ^ = 0, for instance, might have computational advantages one over the other. 
Continuing the integration of V as in the previous method, it is straightforward to check that the 
final answer is given by (69). 

With the experience which we just gained from the example considered above, it is not hard to 
guess that, in general, there is a straight relation between the number of independent zero-modes 
of the Hessian and the number of constraints. Consider the case in which there are M independent 
constraints given by 

^m{A)=0, (70) 

with m = l,2,...,M. We assume that the constraints are expressed such that the gradients of the 
constraints are linearly independent on the constraint surface (for more details, see the regularity 
conditions on how to express the constraints [19]). 
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Since the primary constraints ^m(^) are identically null when expressed as functions of the 
variables, 

Then, for each independent constraint 4>m{'4^i) there is a corresponding independent zero-mode of 
the Hessian given by 

(^)" = ^- P^' 

In other words, upon transforming a Lagrangian given by f(X^), with i = 1,2,...,N, into an 
equivalent one given by i^iX^ — V{^pi), the number of independent auxiliary scalar fields that 
appear in the potential V{ipi) is equal to the rank of the Hessian of f{X^). 

In particular, if / depends on R and E only, the condition of the degeneracy of the Hessian 
matrix reduces to 

fRRfsE - Ire = ■ (73) 

The principal difference between the behaviour of solutions in this special case and in the general 
one has been already noticed when studying small inhomogeneous perturbations on a Friedmann- 
Robertson- Walker (FRW) background [14]. 

A relevant general observation is that we only need the potential V computed on the constraint 
surface, since in general [19] 

V = y|^=o + A0, (74) 

where F|</,=o is the potential on the constraint surface, but the ("primary") constraint (f) identically 
vanishes when expressed as a function of R's. Thus, it does not make any difference whether one 
uses either / = — F or / = — F|()!)=o- 

IV. EXAMPLE: CONDITIONS FOR EXPONENTIAL INFLATION 

As an illustration of the equivalence theorem from the previous section, let us formulate the 
conditions for the existence of a strictly exponential cosmological solution. Consider the action 

Seq = j d^x^{i^iR + ^2^R + ^^C'' + ME-\^R)-V{il^i,il^2,Mi)} (75) 
and the metric 

ds'^ = gi^^dx'^dx'' = a^{ri){dr]^ - df) , a(r?) = e^^*?) , (76) 
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where 

= — + r'^de'^ + sin^ Qd<^^ . (77) 

1 KT 

The conformally transformed metric has the form 

Qi^^ = diag (l, -j^^, -r^sin^^) (78) 

In this section we will restrict our attention to the spatially flat k = case. Then an exponential 
solution produces the de Sitter space-time (it is not so for A; 7^ 0). 
Let us consider the variational derivatives 6Seq/Sijji. 

R - ^nV'l,^2,V'3,V'4)=0 

OR - -^V{i^l,^P2,i^3,^i) = 
dip2 

C - Tr-nV'l,^2,V'3,V'4)=0, 

{E-'^UR) - Ay(^,,^2,V'3,V'4) = 0. (79) 

Furthermore, we need the equation for the metric, which can be obtained by taking the derivative 
of the action with respect to a, 6Seq/6a = 0. This gives 

- 3e-2'^« + 2a'Vi + '^i'W + 2a'^i^i) + 2e-^"Vr - 2y(Vi, ^^2, V'3, ^^4) = 0. (80) 

If we assume that there exists an exponential (in terms of physical time) solution a = aoe^°*, it 
can be inserted into the equations (80) with the following output: 

Furthermore, the equation for the conformal factor has the form 

- 4V - 24i?^Vi - 18i?oV^i - eV'i - 72V^2i^o - 78i/oV'2 

+ 36ifoV'2 + 6^2 + 24iJ^^4 + 44ifo^V'4 + 24ifoV'4 = (82) 

Finding of dS-like solutions in this way implies resolving the system of equations (81), (82), but it 
is not clear whether this can be done in a general form. It should be noted that the r.h.s. of the 
equations (81) - (82) are valid only assuming that we have already used some (unknown) solutions 
for the auxiliary scalars i/j^ there. 
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Let us try another approach for de Sitter solutions. Consider first, as a heat-up exercise, the 
theories with the actions (1) and (2). We know from the Sect. 2 that (1) and (2) are equivalent, 
provided that the functions = f'ji{R) and R = V^{'tp) are inverse functions. Let us use this fact 
to obtain the criterion of dS-like solution for the theories (1) and (2). The dS solution means, in 
the new frame, that 

Rixi^ = ^ Rgixi^ and i?;^i.a/3 = R{9iJ.agvp - Quag 11^) , also R = const . (83) 

How can we see whether the solution of the form (83) is possible or not for the given theory? 
In case of (1) one can take variation with respect to the metric, dgfj,i, = h^i,, and arrive at the 
equation 

^ / - R^.uf'R + (V^V. - g^,ua)f'R = . (84) 

Inserting (83) into the last equation, we arrive at the well known algebraic equation which roots 
R = const give us dS solutions of f{R) gravity: 

lg,.uf-\Rg,.ufR = o =^ RfR = V- (85) 

Now we can do the same for the (2) version of the same theory. The equations equivalent to 

(84) in this case have the form 

]^g^u{i^R-V)-i:R^,^ = Q, also R = V^. (86) 

Inserting (83) into (86), we arrive at the algebraic equation for constant values of at all possible 
dS solutions: 

i^V; = 2V. (87) 

It is fairly easy to see that Eq. (87) is nothing else but the mapping of the final relation in Eq. 

(85) . For this end one has to just use our main relations ^ = fjiiR) and R = V^{ip) in (85). 

One can note that the relation (87) plays exactly the same role for the theory (2) as the relation 
(85) does play for the theory (1). 

The next task is to obtain similar relations for the theories (28) and (29) . One can immediately 
notice that under the conditions (86) we have Xi = BR = and X3 = = 0. Therefore the 
existence of the solution of the form (86) concerns only the dependence of f{Xi) on X2 and X4 
in one case and the dependence of V{'tpi) on ^2 and ^"4 in another case. In all cases we can also 
consider E instead of E. 
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First we deal with Eq. (28). In taking variations of the metric we have to remember that 
after that we shall integrate by parts and then use the conditions (86). Therefore, all covariant 
derivatives, either acting on the variation of the metric h^i, or on curvature tensor components, 
can be safely neglected. In this way we obtain 

SSi = J d^x^{^hf + f'ji-ShR + f'E- ShE] . (88) 

A very simple calculations give 

dhR 

where we denoted h = h^^. Then, for all dS solutions of this theory, the constant invariants R and 



= — - Rh , and 5hE 
dS 4 



= -^R^h = -l-Eh, (89) 
dS 12 2 ' ^ ^ 



E, related by the eonsistency condition E = R^/A in this case, should also satisfy the "on-shell" 
algebraic equation 

f = ^Rf'j, + EfE. (90) 

The last equation is a direct generalization of (85) and has the same theoretical status for the more 
general theory (28) as (85) has for the theory (1). It follows from expressions presented in e.g. 
[14, 20] (though we have been unable to find a paper where it was written explicitly). 

The next step is to obtain the extension of Eq. (87) for the more general case of the theory 
(29). We take all three necessary variations and get 

SS2 „ „ dV 



= ^ R = V( 



SS2 



=^ \{V ^^l^2R-^AE)=--^i:2R-\i^AE, (91) 



dS 

It is easy to see from Eq. (91) that one equation for constant values of 'i/'j, i = 1,2,3,4, at the de 
Sitter solutions (83) has the same form (87) (with ipi instead of ijj and V depending on all ipi) even 
if the Gauss-Bonnet and other terms are present (of course, this does not mean that these terms 
play no role here). 

The other three algebraic "on-shell" equations follow from (81) by excluding Hq (still to be 
found) and using the relation between R and E for a dS solution: 

^^_n -n dv_ii'dvy 
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V. STRING-INSPIRED CASE 

Now let us try to generalize the consideration given above to more general gravitational actions 
containing more complicated terms constructed from the scalar curvature R, the Ricci tensor R^i, 
and the Riemann tensor R^xva^- In general, the corresponding theories can not be reduced to the 
metric-scalar models. For example, the R^^yaP'^^^^ actions involve higher derivatives not only in 
the spin-zero sector but also in the spin-two one. Hence one can expect that the reduction to 
second order equations would require introduction of tensor compensating fields. This is definitely 
true in general, however there is an interesting possibility yet. Let us consider a special space-time 
metric with a restricted number of degrees of freedom, such as the cosmological FRW one. In this 
case we have only one component of the metric - the scale (conformal) factor, which depends on a 
single variable (e.g. conformal time). Then we meet a much simpler situation than in the general 
case, because the tensor structure of the /{R^i^ap) action becomes irrelevant. It might happen 
that the reduction to the metric-scalar theory will be possible in this case. Indeed, this reduction 
concerns only the dynamics of the conformal factor of the metric. After this dynamics is described 
in terms of an appropriate metric-scalar theory, one has to explore other, more complicated aspects 
(e.g. metric perturbations) in the framework of the original higher-derivative theory. However, the 
equivalence with the metric-scalar model may be a useful tool for dealing with a homogeneous and 
isotropic cosmological solution. Therefore, it deserves our attention. 

Let us consider an effective low energy action of (super)string theory (see, e.g. [21]) depending 
only on metric. For the sake of simplicity we assume that the dilaton and effective fermionic degrees 
of freedom (after compactification to o? = 4, because we are interested in the four-dimensional 
physics here) are all in the vacuum state and only metric possesses dynamical behaviour. This 
assumption is perfectly consistent with the known form of the gravitational interaction which 
corresponds to General Relativity (GR), the theory where only massless graviton is propagating. 

In the lowest, first order in the string parameter a', we meet the standard Einstein-Hilbert 
action for gravity. In the next order in a' we meet a set of higher derivative terms, namely 

Rluafs, R% and R\ (93) 

It has been noticed 25 years ago by Zwiebach [22] (see also consequent investigation of the problem 
in [23]) that the choice of the background string metric can be always done in such a way that 
the higher order corrections do not generate unphysical propagating massive ghosts - a typical 
phenomenon for a wide class of higher derivative gravity theories [24, 25]. In the second order in 
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a', the ghosts do not show up if the higher derivative terms (93) enter the following combinations 

E = Rluap-^R% + R\ R". (94) 

Possible metric reparametrizations in the O (ct'^) order have the form 

9nu — > g'nu = 9nu + a' {xi R^,u + X2R g^v) + - , (95) 

where xi^2 s^re arbitrary parameters. The same procedure can be used in higher orders. The above 
transformations can change coefficients (in particular eliminate completely) of all those terms which 
depend on the Ricci tensor or on the scalar curvature i?, only those terms which are constructed 
exclusively from the Riemann tensor may not be modified. The next question is what are the 
physical constraints for those terms which can be modified. As wc have already seen above, the 
E? term can be traded for a scalar field. The first term in the last expression is nothing else but 
the integrand of the Gauss-Bonnet topological invariant (Euler characteristics) of the space-time 
manifold. It is remarkable that this term does not influence the dynamics of the universe. Hence 
one can completely eliminate the relevant 0{R^) corrections to the string effective action by means 
of the metric reparametrization. At the same time, this procedure is not uniquely defined. If we 
require the absence of unphysical ghosts, we may eliminate or not the relevant term depending 
on our own will. The related ambiguity may affect the cosmological solutions [26] and can not be 
fixed without experimental verification (see also corresponding discussion for the case of the gravity 
with dilaton and torsion in [23] and [27]). The R + R^-type action corresponds to the choice of the 
metric parametrization described above. 

In the next order in a', we meet corrections which are cubic in curvature tensor and also the 
i?ni?-type terms. Since the latter term generically leads to the appearance of a ghost scalar [7, 28], 
we investigate the most general metric parametrization in this given order of the former term only. 
The corresponding effective action has the form 

83 = J d^V^ { xiR^ + X2RR,,,R^'' + x^RR^r.apR'''"^^ 

+X4Ri,,R''"R^ + x^R^^apR'^^'^R^" px } • (96) 

Let us now consider the possibility of reduction of the above theory (96) to the metric-scalar 
model in case of the special conformally flat metric g^,^ — a^{r]). Here i] is the conformal time 

and g^i, is the time-independent homogeneous and isotropic metric (78). For our purposes it is 
better to use the variable (j{ri), defined in (76). The transformation rules for the curvature scalar 
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and tensors has the form 



R = e-'^" [R - Q{Vaf - GDa] , 
Rnu = Rf,u-2{Vi,V^.a)-g^^,{na) + 2{Vf,a){VuCr)-2g^4-Vaf , 



(97) 
(98) 



(99) 



For the metric of interest we arrive at the relations for the non-zero components of the above 
curvatures 



R = e 



-2a 



R - 6a' - 6a" 



R, 



-3cr" , Rij =gij[^R- a" - 2a''^ 



1 



(100) 
(101) 

(102) 



where R = —6k = const. 

Using the relations (100) - (102), after some algebra we can rewrite the elements of the action 



(96) in the form 



-gR^ = 



R - 6a'^ - 6a" 



(103) 



gRRauR"" = e 



2a 



R^ - 2R? {2a" + Sa") + 12i? (2a"'' + W + 3a'' J 



72 (a"^ + 2a"' a'^ + 2a" a'^ + a'' 



(104) 



gRR^.^pR'^'"'^ = e 



2a 



^ i?3 - 2R'' (a" + 3(7'') + 12i? (2a"' + 2a" a'' + Sa'^' 



— I Z\a +(j a + a a + a I 



(105) 



gRf^.R'^'R'i = e 



2a 



^R^-R" (a" + 2a") + 3R (a" + 2a")' 

): 



6 ( 5a" V 3a" W 6a"a'' + 4a'' 



'^R^ua^R'^^'^R^'' pA = e^" \\ R^ - 2R?a" + 12i?a'' - 24 (a"' + a") 



(106) 
(107) 
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The sufficient condition of tfie reduction to a metric-scalar tfieory can be easily found by using 
analogy with the ^/^R^ case considered in the previous sections. In the general case this condition 
has the form 

r 1 3 

{xiR^ + X2RR% + X3RR%^^ + X4Rl^ + x^R^^^) = e^" \yiR^ + y2a" + j/a^^'^J , (108) 

where j/1,2,3 are some additional arbitrary coefficients. Thus we obtain 10 algebraic equations for 
the 8 variables xi^2,3,4,5 and 1/1,2, 3- In fact, the number of the variables can be immediately reduced 
to 7 by noticing that the coefficient X5 can not be made zero in the string induced gravity by 
means of the metric reparametrization, while all other coefficients can [22]. Therefore without 
losing generality we can set = 1. Now, as far as the number of equations is much greater than 
the number of independent variables, it is not certain that the solution of these equations exist. 
As we shall immediately see, it exist only for some particular, but the most relevant case. Let us 
remember that our prime interest is the inflationary epoch, where the space curvature is negligible. 
Then, as a first step we can look for the solution in the simplest A; = case, where ^ = and the 
number of equations is even smaller than the number of independent variables. In this case the 
equations become 

-yi = 216x1 + 72x2 + 72x3 + 30x4 + 24, (109) 
-32/iy3 = 648x1 + 144x2 + 72x3 + 18x4 , (110) 
-3y22/i = 648x1 + 144x2 + 72x3 + 86x4 , (111) 

-y| = 216x1 + 72x2 + 72x3 + 24x4 - 24. (112) 

Using the pairs of equations (110), (111) and (109), (112) we obtain the relation 

2/32/2(2/3 - 2/2) = (ys - 2/2)(2/3 + 2/3^2 + 2/2) = 6x4 . (113) 

These equations can be satisfied only for 2/2 = 2/3 = 2/ and X4 = 0. Let us remark that the constraint 
2/1 = 2/2 means that possible dependence on the conformal factor for the theories reducible to metric- 
scalar models can be only the same as for the ^/^R^-c^lSe. All other choices are not reducible. 
The general solution corresponds to 
X 2 2/*^ 

2^1 = -3 - g + 2jg , a;2 = -2x-l + — , X3 = x , X5 = 1 , (114) 

where x,y are arbitrary parameters. It is easy to see that if the constrains (114) are satisfied, we 
have the same reduction to a metric-scalar model as in the case of the higher derivative term 

,,3 



y 

216 



Jd'^x^/^R^. (115) 
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At the level of the metric-scalar model we have, in the string-induced case, exactly the same 
situation as in the case of the term (115). 

For the sake of completeness, we consider the more complicated A; 7^ case. An important 
observation is that the constraints (114) must hold also for k 0. Then elementary analysis shows 
that in this case there are no solutions of eq. (108). Therefore, mathematically the reduction 
to the metric-scalar model can not be exact in the general case, but only an approximate one. 
However the quality of this approximation is indeed excellent because the role of k during inflation 
is negligible. 



As another illustration of the effectiveness of our approach, consider the Eddington-like action 
of gravity [29, 30] [48], 



Here and in what follows we use the notation \Rnu\ = \ det (i?^i/)|, furthermore a is a dimensionless 
parameter and R^i, is the (symmetric) Ricci tensor constructed from the symmetric affine connec- 
tion r^^. The action (116) is equivalent to the vacuum Einstein-Hilbert action with a non-zero 
cosmological constant [31, 33]. In what follows wc will show it is straightforward to achieve the 
same result by using the method we have considered in the previous sections. Moreover, in the 
known approaches [33, 34], it is common to consider the first order formalism for gravity, taking 
r^j^ to be independent from the metric g^^,. Our method does not require this restriction and, 
moreover, R/^i, can be traded for any other symmetric tensor, e.g. to some combination of the 
torsion fields (e.g., the one considered in [27]). Of course, the corresponding model will not be 
equivalent to the GR, but it can be mapped to a dual theory in the same way as we will describe 
below. 

The dual equivalent action should have the form 



where J'*^ is an auxiliary field. One can immediately note that there is a unique functional form 

of the potential function V{ J) = V{ J^'^) which is compatible with the covariance of the action 
(117). Indeed, the covariance requires that J'^'^ should be a tensor density and there must be such 
symmetric tensor quantity ^^u, such that 



VI. METRIC DUAL FOR THE EDDINGTON-LIKE GRAVITY 




(116) 




(117) 



(118) 
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where • ^ai^ = and $ = | det This relation (118) can be easily inverted, so we get 

$ = det = , } , and ^f"" = , ^ J'^'^ . (119) 

^ det(JM-) 71 det (J'^'^) I ^ ^ 

As a result, we arrive at V{Jf^'^) = k ■ \/$, where k is some constant. It is easy to see from the 
relations (119) that this means 

y( j^^) = ^ ^ . (120) 

Now let us see whether we can arrive at the same result (117), (120) by using the method 
described in the previous sections. By taking derivative of the function f{R^v) = o^^/\R|^u\ and 
inverting it, we obtain 

The equivalence of the two formulations requires that J^^ = J^^{Rap), R^v = Rij,u{J"^) and, 
also, that 



J'^'' ■R^,-ViJn = a^\R,^A- (122) 

Now, taking the partial derivatives d/dJ^^ of the equality (122), after some simple algebra we 
arrive at 

j'''=%^\{K'r ^ (123) 

where {^R~^)'^^ means the matrix inverse to R^i,. It is easy to obtain the relation 

(R-^f- ^ « J^-^ (124) 

which shows that J''^ is a tensor density. The corresponding tensor field is defined through the 
relation (118). After some small algebra we arrive at the equation 

Finally, integrating this equation we get the expected result 

F(J'*'^) = , ^ = - (126) 

that is nothing else but Eq. (120) with k = A/a. 

It is clear that the consideration presented above opens the way for some interesting applications. 
One can, at the first place, identify the auxiliary tensor quantity directly with the space-time 
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metric g^jj, but this is not the only one possible choice. Let us note that the real identification of 
the metric occurs when matter is introduced into the theory. One can assume, for instance, that 
there is a scaling relation = A • g^i, or even more complicated one, like = B{x) ■ g^i,, where 
the scalar field B{x) depends on the space-time coordinates. The corresponding generalization of 
the Eddington-like gravity theory looks interesting and perhaps deserves further discussion. The 
method presented here can be immediately extended to the case when R^^, in the initial action 
(116) is traded for some other tensor, e.g. for + Cg^u- In this case the auxiliary field ^^i, is 
naturally identified with the second metric and we arrive at the bimetric theory of gravity [31]. 

VII. A BRIEF NOTE ON QUANTUM (NON) EQUIVALENCE 

In the previous sections 2-4 and 5 we have presented several examples of classically equivalent 
theories. It looks interesting to see what happens with this equivalence at the quantum level. 
It is well known that the quantum equivalence does not imply the quantum one. One can find 
statements of this type in Rcfs. [9, 11] and also in [17]. Of course, classical equivalence leads to 
quantum equivalence at the level of tree diagrams or imaginary parts of one-loop diagrams which, in 
particular, describe creation of real particles and field perturbations by external gravitational fields 
and do not need renormalization. That is why primordial spectra of scalar and tensor perturbations 
generated during inflation in f{R) gravity (1) (in the models which admit it, say, in the R + B? 
model [35]) coincide with those for inflationary models with a non-minimally coupled scalar fleld 
in the limit of a large negative coupling |^| S> 1, see e.g. [36-38], since the latter models reduce to 
the form (2) after neglecting the kinetic term of the scalar field that is justified in this limit. Note 
that the Higgs inflation [39] belongs to this class, too. 

The problems of establishing general quantum equivalence become much more complicated 
when taking into account the need for renormalization. For the sake of simplicity we consider the 
theories (1) and (2), but the generalizations to other cases are indeed possible. 

The quantum equivalence at the one-loop level means, at the flrst place, that the relations 
like V' = f'{R) E^nd R = V'{il^) do hold for the one-loop counterterms. In reality, this is not 
necessary so, because these relation may also require renormalization. Moreover, even if the one- 
loop divergences do satisfy this requirement, it is very unlikely that some sort of relation between 
the two (classically equivalent) theories will hold beyond one-loop approximation. 

The main difficulty of discussing the equivalence of the two formulations, e.g., (1) and (2), or 
(27) and (28) - is the fact that the corresponding theories are not renormalizable. At the same time. 
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even one- loop divergences may be non-equivalent, as it happens with the tensor-scalar model in 
different conformal frames [40, 41]. Indeed, the non-equivalence here means that one can not easily 
find an explicit transformation which would link the two expressions for one-loop counterterms 
in the two representations. At the same time, since both theories are non-renormalizable, one 
can speak about certain qualitative equivalence in a sense that, in both frames, the necessary 
counterterms have the structure distinct from the one of the initial action. 

One can look from another side and compare the UV completion of the two, classically equiv- 
alent, theories. It is easier to perform such discussion for a more general cases of the theories 
(27) and (28). The theory (28) with "frozen" scalars ■0* is renormalizable [25] (see also [42] for a 
more detailed introduction). Therefore, in order to construct renormalizable theory out of (28), 
one has to complete the action by certain second- and fourth-derivative terms constructed from 
the scalars ip'^. In case of one scalar field this completion has been considered in [43] and more 
recently in [44] in relation to inflation. The generalization to the many-scalar case is obvious, but 
the output would be quite cumbersome. The situation with the general higher-derivative model 

(27) is much more complicated. First of all, there is a very strong difference between polynomial 
and non-polynomial functions f{Xi). In the last case the problem of quantum formulation is un- 
clear (despite potentially interesting, see e.g., [45]). Contrary to that, in the polynomial case the 
prescription for constructing renormalizable and super-renormalizable theories of quantum gravity 
is known [24]. If the highest power of curvature tensor in f{Xi) is iV > 3, one has to introduce 
into Lagrangian all possible covariant term of this dimension, including the additional terms of the 
following form: 

7V-2 

J2 {cxkRf.ua'^R'"' + PkRo'^R) ■ (127) 

fe=0 

If the largest order coefficients aN-2 and Pn-2 are nonzero, the theory is super-renormalizable 
[49]. 

Now we are in a position to compare the renormalizability properties of the theories (27) and 

(28) . The UV completions described above correspond to the possible counterterms and, in these 
two cases, these counterterms are dramatically different. For example, the UV completion in the 
case of scalar-tensor theory (28) describes the propagation of only two spin-2 states, namely of 
graviton and of the massive spin-2 ghost. At the same time the UV completion of the theory 
(27) has (for AT > 3) at least one more spin-2 massive particle [24]. Thus, we can conclude that 
the renormalization properties of the two classically equivalent non-renormalizable theories are, in 
general, quite different. In particular, they have very distinct UV completions and, consequently. 
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very distinct structure of couterterms. 

VIII. CONCLUSION 

We have described in the previous sections how to perform a mapping of a gravitational theory 
of the type f{X^) into a theory with auxiUary scalar fields. The number of these fields is determined 
by the rank of the Hessian matrix of /(X*). The scheme which we have described is more general 
than the ones known before. In particular, it enables one to deal with the constrained case and 
leads to an auxiliary-field representations for the actions which look like renormalization-group 
corrected vacuum actions in gravity theories like Eq. (49). On the top of this we have formulated 
the general conditions for the exponential expansion of the Universe and extended the analysis for 
the string-inspired case, where the treatment with auxiliary scalars is possible only for a FRW-like 
solution. Finally, we have shown that our approach is perfectly applicable to the Eddington-like 
gravity models and to a wide class of their extensions. Depending on the initial model, our approach 
shows cither how the metric emerges in a theory which had, initially, only affine connection, or 
leads to a bimetric theories of gravity. 
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